function [data, data_format] = readfeaturefile(filename)

fid = fopen(filename,'r');
fgetl(fid);
fgetl(fid);

data_format = fgetl(fid);
C  = textscan(fid,'%f');
data = C{1};
data = reshape(data,15,[])';

% format: x y z node_id node_label [features]
fclose(fid);